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Abstract. The influence of small random perturbations on a deterministic dynam¬ 
ical system with a locally stable equilibrium is considered. The perturbed system is 
described by the ltd stochastic differential equation. It is assumed that the noise does 
not vanish at the equilibrium. In this case the trajectories of the stochastic system may 
leave any bounded domain with probability one. We analyze the stochastic stability of 
the equilibrium under a persistent perturbation by white noise on an asymptotically 
long time interval 0 < t where 0 < /i <C 1 is a perturbation parameter. 
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1. Problem statement 


Consider the system of ordinary differential equations: 


( 1 ) 


(j'V 

— = f(x,f), t>0, xG 


Let X = 0 be an equilibrium, f(0,t) = 0. Suppose the vector-valued function f(x, t) = 
(/i(x, t),..., t)) is continuous, satisfies a Lipschitz condition: |f(x, t) — f(z, t)\ < L|x — z| 

and a growth condition: |f(x, f)| < M(1 -|- |x|) for all x, z G M”, f > 0 with positive constants 
L, M. Assume that there exists a local Lyapunov function U(x, t) G x M) satisfying 

the following properties: 


( 2 ) 


dt 


< U(x,f) < A|x|^ l^xUl^ < R|x|^ <C, 

defdV ^dV w W| I . n 

+ <-W V|x|<ro, t > 0, 


( 1 ) 


dt 


k=l 


dxi 


with constants A, B,C,''y,rQ > 0. This implies that the solution x(t) = 0 is locally asymptot¬ 
ically stable (see [1]). In this case, system (1) may have other stable fixed points. Note that 
such Lyapunov functions are constructed in the study of nonlinear differential equations (see, 
for instance, [2,3]). In this paper we investigate the stability of the trivial solution under a 
persistent perturbation by white noise. 

Together with (1) we consider the ltd stochastic differential equation: 

(3) dy(t) = f(y,t)dt-FpG(y,t)dw(t), y(0) = yo G M”. 


Here w(t) = (tci(t),..., Wmit)) is m-dimensional Wiener process defined on a probability space 
(r2,J^, P), G{y,t) = {dijiyA)} is a given continuous n x m matrix, satisfies the Lipschitz 
condition: 


(4) ||G(y,t)-G(z,t)|| < L|y-z| Vy,zGM", t > 0, 
and the growth condition: 

(5) ||G(y,t)|| =^max| 5 fy(y,t)| < M(1-h |y|) VyGM'", t > 0. 
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A small positive parameter 0 < p <C 1 controls the intensity of the random perturbation. 
The matrix G{y,t) and the initial data yo are assumed to be deterministic. Note that these 
constraints are sufficient for the existence and uniqueness of a global solution to the initial 
value problem (3) for all yo G M” (see [4, §5.2], [5, §3.3], [ 6 , §6.2]). In addition, suppose that 
the trivial solution y{t) = 0 is not preserved in the perturbed system, namely, G(0, t) ^ 0. The 
goal of this paper is to describe a class of matrices G guaranteing the stability in probability 
of the trivial solution to system ( 1 ) under a persistent perturbation by white noise. 

Remind the concept of stability in probability [5, §5.3], [7] or stochastic stability [ 8 , Ch. 
2], [9, Ch. 11]: 

Definition 1. The solution x(t) = 0 to system (1) is said to be strongly stable in probability 
with respect to white noise uniformly for G from the set V, if 

We,n>0 3(5, A>0: V|yo|<(5, /i < A, GeV 
the solution y{t) to the initial value problem (3) satisfies the inequality: 

P(sup |y(f)| >e)<i'. 

t>o 

Thus, if the trivial solution is strongly stable, then the trajectory of the process y(t) with 
sufficiently small initial values and perturbations does not leave an arbitrary small neighbour¬ 
hood of zero with probability tending to one. If the solution x(t) = 0 is weakly stable^, then 
almost all realizations may leave the ball = {y E M"" : |y| < e}, but at each instant t = t^, 
the random variable y(f*) satishes |y(t*)| < e with probability tending to one. It is easy to see 
that strong stability implies weak stability, but the converse is not true in general [5, §6.11]. 

Many authors (see [5-10]) have considered the stability of stochastic differential equations 
with G(0,f) = 0. However, the problem of stability under persistent perturbations was investi¬ 
gated only in several works [5,11-14]. Let us briefly recall the main known results concerning 
the case G(0, t) ^ 0. Firstly, from [11] it follows that there is no strong stability in autonomous 
systems (where both f and G do not explicitly depend on t): the trajectories of perturbed 
system may leave any bounded domain with probability one. In [5, §7.4], Khasminskii proved 
that if ||G(y,t)|| decays sufficiently fast at inhnity t ^ oo and there exists a suitable global 
Lyapunov function, the solution x(f) = 0 remains stable in a certain sense. On the other 
hand, the existence of a global Lyapunov function for stochastic equation (3) is sufficient for 
weak stability w.r.t. white noise with uniformly bounded matrix G (see [12]). However, if the 
deterministic system has at least one more stable hxed point x* 7 ^ 0 , the solution x(t) = 0 is 
not even weakly stable w.r.t. such class of perturbations. Thus, the stability under persistent 
perturbation by white noise occurs in a fairly narrow class of systems, and it is therefore worth 
while to consider the problem of stability on a hnite time interval (see [11,13]). One variant 
of this approach is to hnd the largest possible time interval [0; T] on which solutions to the 
perturbed equation are close to the equilibrium of the deterministic system. It was shown 
in [13, Ch. 9] that if the matrix G is uniformly bounded and there exists a local Lyapunov 
function having properties ( 2 ) with 7 = 0 , then the equilibrium is strongly stable on the interval 
0 < f < 0{pr‘^). Similar results on weak stability were obtained in [14], where a barrier for 
the Kolmogorov equation (see [9, §4.4]) was constructed using a Lyapunov function known in 
a neighborhood of the equilibrium. However, the stability of solution on longer time intervals 
remains an open question. The present paper is devoted to solving this problem. 

In this work we prove that if G is uniformly bounded matrix, then the locally asymptotically 
stable solution x(f) = 0 to the deterministic system is strongly stable with respect to white 
noise on an asymptotically long time interval 0 < f for any natural number iV > 1. 


^In this case it is assumed that P(|y(t)| > e) <v for all t > 0. 
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2. Stability on asymptotically long time interval 
Let us clarify the definition of stability we shall use in this section. 

Definition 2. The solution x(f) = 0 to system (1) is said to be strongly stable in probability 
with respect to white noise on a time interval [0; T] uniformly for G from the set V, if 

Ve, i^>0 3(5, A >0: V|yo|<(5, /i < A, GeV 

the solution y{t) to the initial value problem (3) satisfies the inequality: 

P( sup |y(t)| >e)<u. 

te[0;r] 

We dehne the class ^ as a set of matrices G{y,t) satisfying estimates (4), (5), and 

II ^ +MI ^ def G ■ G* 

sup ||a(y,t)|| < oo, a = —-— = 

|y|<ro,t>0 ^ 

Let Ah be a set of matrices G from the set A such that ||(j(y,t)|| < h for all |y| < tq, t > 0. 
Let us also dehne the operator (see [5, §3.6]): 

n n 

C dt + '^ fk{y, t)d^^^ + 

k=l *,J=1 


We have 


Theorem 1. Suppose that unperturbed system (1) has a Lyapunov function ld(x, t) satisfying 
(2). Then for all N E N, h > 0, and 0 < x < 1 the solution x(t) = 0 to system (1) is strongly 
stable in probability with respect to white noise on the time interval [ 0 ; uniformly for 


G E Ah- 


Proof. It is reduced to constructing a suitable Lyapunov function for stochastic differential 
equations (3) on the basis of the Lyapunov function V (x, t) for unperturbed deterministic 
system ( 1 ). 

Let h > 0, 0 < X < 1, > 0, and 0 < e < tq be arbitrary positive constants, y{t) a solution 

to system (3) with |yo| < S, G E Ah, and tx a hrst exit time of the solution y{t) from the 
domain 

X {(y,t) G : |y| < e, 0 < t < T}, T > 0. 

We dehne the function st = min{ri, t}. Then y(st) is the process stopped at the hrst exit time 
from the domain X. 

Let us hrst consider the case N = 1. Here we reproduce the arguments of [13, Ch. 9] with 
some modihcations. The Lyapunov function for system (3) is constructed in the form; 

( 6 ) Vi{y, t-,T) = V{y,t) + p^n^hG ■ {T -t), 


where V{y,t) is the Lyapunov function for system (1). It is easy to see that Vi>V>0 and 


CVi = 


dt 


+ /i^ ^ aij d^^d^.V - p^n^hG < -yH < 0 


for all (y,t) G X. Hence Vi(y(st),st) is a nonnegative supermartingale [5, §5.2]. From the 
properties of the Lyapunov function V and the dehnition of the function St, we get the following 
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estimates: 


( 7 ) 


P( sup |y(t)| > e) 

0<t<T 


< 

< 


P(sup |y(t)r>£^)< 

0<t<T 

P( sup V{y{t),t)>e^)< 

0<t<T 

P( sup Vi{y{t),t]T) > = 

0<t<T 

P{sup Vi{y{st), St; T) > < 

t>o 


Vi(yo,0;T) 


The last estimate follows from Doob’s supermartingale inequality [15, §7.3]. Note that the 
initial point yo is assumed to be deterministic. We dehne T = and 


<5 = 


\2AJ 


1/2 


^ \2n^hc} ■ 


Then Vi(yo,0;T) = A|yop + < e^v for all |yo| < 6, p < A. Taking into account (7), 

we obtain the estimate: 


(8) P( sup |y(t)| >e)<u. 

0<t<T 

This implies that the solution x(t) = 0 is strongly stable as t 

Now let us prove that the trivial solution of system (1) is stable as t < . The Lyapunov 

function is constructed on the basis of the functions V and Vp. 

V 2 {y, t; T) = (P(y, t)f + p'^aMy, t; T), 

where parameters Oi and T are defined later. The action of the operator C on V 2 satishes the 
inequality: 


< -2yV^ + [2 n^h{B + C)- uiq] V 

as |y| < tq. We choose Oi = 2 n?h{B + C)'y~^, then V 2 >V‘^>Q and CV 2 < 0 for all (y, t) G X. 
Hence V 2 {y{st), St;T) is a nonnegative snpermartingale, and the following inequalities hold: 


(9) 


P( sup |y(t)| > e) 

0<t<T 


< 

< 


P(sup |y(t)r>£^)< 

0<t<T 

P(sup {V(,y(t),t)f>E*)< 

0<t<T 

P( sup V2{y{t),f,T) > e^) = 

0<t<T 

P{sup V 2 {y{st), St; T) > e^) < 
t>o 


V2{yo,0;T) 


At this step we choose T 


jj, and define the parameters <5 > 0, A > 0 such that 


3^254 

2 



+ C < e^v. 


Then the inequality V 2 (yo 5 0;T) < holds for all |yo| < 6, fi < A and G G Ah- Combining 
this with (9), we get estimate (8) showing the stability on the time interval [0; 

Finally let us prove the stability on the interval 0 < t < with N > 3. In this case 

we construct the perturbed Lyapunov function in the following form: 


V^{y,t;T) = {V{y,t)f + fiW-iVN-i{y,t;T), 

W(y,f;T) = {V{y,t)T + fi"ak-iVk-t{y,t-,T), k = 2,...,N-l, 







WHITE NOISE PERTURBATION OF LOCALLY STABLE DYNAMICAL SYSTEMS 


5 


where the function i/i(y,t;T) is determined by formula ( 6 ) and T = The additional 

condition £14 < 0 , imposed on each function 14 in the ball |y| < vq, leads to the choice of the 
parameters a^. Let = (fc + l)n^h{B + C)'y~^, then the following inequalities hold: 

£14 < - 371 ^^ + 2/^2 (sn^h (B + C)- 027 ) < 0, 


£14v < -N-fV^ + {N- l)/i2 (^Nn^h (B + C)- aiv- 17 ) 1^^“^ < 0 


for all |y| < tq, t > 0, and = 3,4,.... Hence the function 14v(yo('St), St] T) is a nonnegative 
supermartingale. Since Idv > > 0 for all (y, t) G X, we have the inequalities: 


( 10 ) 


P( sup |y(l)| > e) 

0<t<T 


= P(sup 

0<t<T 

< P(sup >e^^)< 

0<t<T 

< P( sup l/Ar(y(t),t;T) > = 

0<t<T 


= P(supl4r(y(st),St;T) > £ )< 

t>0 


2N\ / 14r(yo,0;T) 


r2N 


It can easily be checked that l/Ar(yo,0;T) = O(|yop'^) + as |yo| —)■ 0 and /i —)■ 0. This 

yields that there exist <5 > 0 and A > 0 such that 14v(yo) 0 ; T) < for all |yo| < 5, /i < A. 
Combining this with (10), we get estimate ( 8 ). Therefore, V A G N, h > 0 the solution x(f) = 0 
to deterministic system ( 1 ) is strongly stable with respect to white noise on the time interval 
[0; uniformly for G G Ah- This completes the proof. □ 


Remark 1. The theorem holds true if we choose T = /i~^^A(|yo|) with a positive contin¬ 
ues function A( 2 :) such that X{z) ^ 0 as z ^ 0. In this case we get the inequality: 0 < 
14v(yo,0;T) < m(|yop^ + + A(|yo|)) as |yo| < ro and /r < 1 with a positive constant m. 

The parameters 5 and A are chosen such that 6“^^ + A((5) < e^^v/{2m), A^^ < e^^i'/[2m). 
Hence we have (8) for all \yo\ <6 and fv < A. This implies that the solution x(f) = 0 is stable 
on the asymptotically long time interval 0 <t < yU“^^A(|yo|). 


Remark 2. If the Lyapunov function P(x, f) have properties (2) with 7 = 0, then the proposed 
construction of the perturbed Lyapunov function guarantees the stability of the trivial solution 
only on the time interval 0 <t < 


3. Stability under damped perturbations 


Let us consider a narrower class J' consisting of matrices G{y,t) satisfying (4), (5), and 
having at least one continues function ((t) G Xi(M+) such that 

sup ||(T(y,t)|| < C(t) Vf>0, a ^ . 

|y|<»’o ^ 

Let Jh be a set of matrices G ^ J such that 

00 

j C{0) dO < h. 

0 


We have 


Theorem 2. Suppose that unperturbed system (1) has a Lyapunov function l/(x, f) satisfying 
(2) with 7 = 0 . Then for all h > 0 the solution x(f) = 0 to system (1) is strongly stable in 
probability with respect to white noise as t >0 uniformly for G ^ Jh- 
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Proof. We fix parameters /i > 0, z/ > 0, and 0 < £ < rg. Suppose y(t) is a solution to 
system (3) with |yo| < <5, G G jTh, Te is a first exit time of the solution y(t) from the domain 
Be = {y eR^ : |y| < e}, and s* = min{r£,t}. Then the function y(st) is the process stopped 
at the hrst exit time from the ball 

The Lyapunov function for system (3) is constructed in the form 


Vfiiy, t) = v (y, t) + pVG • [h - 



From the dehnition of the set Jh and properties of the function V it follows that > 1^ > 0 
and 


^ dt 


( 1 ) 


+ ~ n^'nfCCit) <0 V |y| < ro, t> 0. 


* 3=1 


Consequently, the function V^(y(st),St) is a nonnegative supermartingale, and we have the 
estimates: 

P(sup |y(f)| > £) = P(sup|y(si)|^ > < 


t>0 


t>0 


< P(sup V (y(st), St) > e^) < 


( 11 ) 


t>o 


< P(sup Vf,{y{st), St)>s) = 

t>o 


= P(sup V';.(y(si), St) > e ) < 


2^ / v;.(yo,o) 


t>0 


The last estimate follows from Doob’s supermartingale inequality. We choose 5 > 0, A > 0 
such that A|yop + jjf'ofCh < for all |yo| < <5, /i < A. Hence, V^{yQ, 0) < Combining 
this with (11), we get estimate (8). Therefore, the solution x(t) = 0 to system (1) is strongly 
stable with respect to white noise as f > 0 uniformly for G E fTh- Cl 


4. Examples 

1. Let us consider the perturbed ltd stochastic differential equation: dy(t) = ydw(t). The 
solution x(t) = 0 to the unperturbed system x = 0 is stable (but not asymptotically stable); 
the Lyapunov function V = satishes estimates (2) with 7 = 0 . Let us show that the trivial 
solution is strongly stable w.r.t. white noise (with G = 1) on the time interval 0 < f < o(/i“^). 

It can easily be checked that the function y{t) = yo + fiw(t) is the solution to the perturbed 
equation. For simplicity, we assume |/o = 0. Using the known formula for the Wiener process 
(see [16, p. 173]), we get 


(2fe+l)c 



c, T = const > 0. This implies that for all 0 < x < 1, e > 0, and u > 0 there exists 
A = (£^/|41n(i//\/8)|)^A such that 

P( sup \y{t)\>e)<u 

as p < A. In this example the strong stability is not preserved as t ^ /i“^. The stretching of 
the time interval can be achieved by imposing the additional constraints on the perturbation 
(see theorem 2). 
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2. Let us consider a more complicated equation: 

dyit) = - dt + dw(t), 0 </i 1 . 

1 + r 

The corresponding deterministic system 

dx x{l — x^) 
dt 1 + 

has three equilibria: {—1, 0,1}. It is easily shown that the trivial solution x(t) = 0 is locally 
exponentially stable; the Lyapunov function V{x) = x"^ satishes the inequality dV/dt < —V for 
all |a:| < l/\/3, f > 0. From [12] it follows that in this case there is no stability of the trivial 
solution under white noise on a semi-inhnite time interval with G & A. However, it follows 
from theorem 1 that the solution x{t) = 0 is strongly stable on the hnite but asymptotically 
long time interval 0 < t < 
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